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The  Madelung  energies  of  simple  metal  spheres  are  calculated  for  bcc,  fee,  and  hep 
structures  as  well  as  for  fully  relaxed  structures.  It  is  found  that  for  clusters  of  up  to 
40  atoms  the  structure  generally  does  not  have  the  symmetry  of  any  simple  lattice.  The 
variation  of  the  structural  part  of  the  total  energy  is  shown  to  be  slightly  smaller  than 
the  variation  arising  from  the  filling  of  discrete  single  electron  energy  levels. 
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'Recent  experiments7^  on  the  ’magic  numbers’  of 
small  sodium  clusters  suggest  that  variation  of  the  binding 
energy  of  clusters  as  a  function  of  their  size  is  dominated 
by  the  filling  of  energy  levels  of  electrons  moving  in  a 
spherical  potential  .  This  simple  interpretation  is  sup¬ 
ported  also  by  good  agreement  of  the  calculated  elec¬ 
tronic  polarizabilities4  with  the  experimentally  measured 
ones5.  Ab  initio  calculations6,7  of  the  structures  of  the 
clusters  have  been  performed  only  for  the  smallest  clus¬ 
ters.  These  give  the  striking  result  that  for  up  to  six  atoms 
the  ground-state  structure  of  the  cluster  is  planar.  How¬ 
ever,  these  nonspherical  structures  are  not  entirely  in  con¬ 
tradiction  with  the  spherical  jellium  picture,  since  the 
structures  can  be  well  understood  by  studying  the  single 
electron  wave  functions  in  a  spherical  potential:  in  the 
case  of  a  partly  filled  shell  the  charge  density  will  be 
nonspherical  and  the  total  electrostatic  energy  can  be 
reduced  by  relaxing  the  ions  to  a  configuration  which  has 

the  shape  of  the  electron  density  distribution. 

\ 

^JFor  large  clusters  it  is  expected  that  the  electron  den¬ 
sity  distribution  starts  to  be  more  and  more  spherical  as 
the  cluster  increases.  It  is  then  appealing  to  make  a 
spherical  approximation  in  determining  the  cluster  struc¬ 
tures.  Recently,  Iniguez  et  al8  have  used  this  approxima¬ 
tion  and  the  density  functional  scheme  to  calculate  the 
binding  energies  of  sodium  clusters  in  some  bcc  and  fee 
configurations.  A  related  calculation  was  published  ear¬ 
lier  by  Buttet  et  al9.  The  results  show  large  variations  in 
the  cluster  binding  energy  as  a  function  of  the  cluster 
size,  and  Inigues  et  al8  conclude  that  the  actual  clusters 
have  "disordered"  structures  which  are  more  spherical 
than  those  that  can  be  formed  from  simple  bcc  and  fee 
lattices. 

■the  purpose  of  this  communication  is  to  show  that  if 
the  electron  density  distribution  is  assumed  to  be  spheri¬ 
cally  symmetric,  the  resulting  ionic  configuration  is  deter¬ 
mined  mainly  due  a  classical  Madelung  energy.  Further¬ 
more  it  is  shown  that  the  ground  state  structures  of  the 
clusters  are  generally  not  small  pieces  of  simple  lattices 
(bcc,  fee,  or  hep),  but  take  always  a  more  spherical  shape 
(e.g.  icosahedral  in  the  case  of  a  13  atom  cluster).  The 
variation  of  the  structural  energy  as  a  function  of  the 
cluster  size  is  found  to  be  somewhat  smaller  than  the 


variation  arising  from  the  filling  of  single  electron  energy 
shells  giving  support  to  the  idea  that  the  dominant  magic 
numbers  are  determined  by  the  electronic  shell  structure 
of  a  square  well  potentiaP^^ 

In  the  density  functioned,  method,  the  total  energy  of 
the  cluster  is  written  as8 

Et  =  7T«1  + 

+  ..hr  + 1 (i) 
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where  T  and  Exc  are  the  kinetic  and  exchange-correlation 
energies,  n(r)  is  the  electron  density,  R,  the  ion  site,  and 
Vps  the  electron-ion  pseudopotential.  Now  we  define  a 
fictitious  charge  density  distribution  n0  as  a  step  function 

"oM  =  n0Q(Rc-r),  (2) 

where  Rc  is  the  cluster  radius,  and  nQ  is  the  average  con¬ 
duction  electron  density  in  an  infinite  metal.  n0(r)  is  equal 
to  the  compensating  positive  background  density  in  a  jel- 
lium  cluster.  By  adding  and  subtracting  electrostatic 
energy  terms  in  Eq.  (1)  we  can  write  the  total  energy  as 

Er  =  T[n]  +  Exc[n)  +  Es( 

1  ri(rMr')-n0(r)/!0(r') 


+  X^/rl^Cr-^Mr)  +  Xj^r— 
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where  the  Madelung  energy  EAi  is  defined  as 
1  -  ,  n0(r)nQ(r') 


The  advantage  in  writing  the  total  energy  in  this  form  is 
that  most  of  the  structural  dependence  of  die  energy  is 
included  in  the  Madelung  energy.  This  is  analogous  to  the 
classical  cleavage  energy  in  calculating  the  surface  energy 
in  the  same  method10.  The  two  last  terms  in  Eq.  (3)  also 
depend  explicitly  on  the  structure.  Using  for  example  a 
Heine-Abarcnkov1 1  or  Ashcroft1-  type  pscudopotential 
which  equals  the  1/r  potential  beyond  a  core  radius,  say 
r  ,  it  is  easy  the  see  that  die  structure  dependent  parts  of 


a 


the  two  last  terms  in  Eq.  (3)  roughly  cancel  each  other. 
These  two  terms  can  be  written  as 


Xj^rV^(r-R,)[«(r)-*0(r)] 
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The  structural  dependence  of  the  first  term  comes  only 
due  to  the  difference  n-n0  which  is  markedly  nonzero 
only  in  a  narrow  region  at  the  surface.  The  structural 
dependence  of  the  second  term  is  also  small  since  the 
integral  is  limited  only  over  the  core  regions  of  the  pseu¬ 
dopotentials.  Since  the  other  terms  of  Eq.  (3)  also 
depend  only  weakly  on  the  structure13  it  is  reasonable  to 
assume  that  only  the  Madelung  energy  depends  strongly 
on  the  structure.  Thus  in  determining  the  structure  of  the 
clusters,  the  calculation  of  the  electronic  structure  is  not 
needed  at  all  (note  that  all  this  is  based  on  the  assumption 
that  the  resulting  electron  density  distribution  as  a  result 
of  a  self-consistent  density  functional  calculation  is  spher¬ 
ical).  This  fact,  that  the  structural  dependence  of  the  total 
energy  is  dominated  by  the  Madelung  energy,  has  been 
noticed  earlier  in  calculating  surface  energies14-16  and 
void  formation  energies16  in  simple  metals. 

We  have  calculated  the  Madelung  energies  of  clusters 
up  to  40  atoms.  The  cluster  structures  considered  are  bcc, 
fee,  and  hep  lattice  structures  and  the  fully  relaxed  struc¬ 
ture  which  minimizes  the  Madelung  energy.  The  results 
for  the  binding  energy  (i.e.  negative  of  the  Madelung 
energy)  is  shown  in  Fig.  1.  The  relaxed  ground  state 
always  has  a  much  higher  binding  energy  than  any  simple 
lattice  configuration.  In  calculating  the  Madelung  ener¬ 
gies  of  bcc,  fee,  and  hep  clusters  the  lattice  constant  was 
fixed  to  correspond  to  the  electron  density  n0.  The  calcu¬ 
lations  were  repeated  for  all  three  lattices  by  taking  the 
center  of  the  cluster  to  be  a  lattice  site,  tetrahedral  inter¬ 
stitial  site,  octahedral  interstitial  site,  midpoint  of  two 
nearest  atoms,  or  the  center  of  a  equilateral  triangle  of 
three  closest  atoms.  The  optimal  lattice  structures  (bcc, 
fee,  or  hep)  for  clusters  from  l  to  40  atoms  are  given  in 
Table  1.  These  structures  are  generally  very  different 
from  the  actual  ground  state  structures,  which  were 
obtained  by  letting  the  ions  relax  to  positions  which 


minimize  the  Madelung  energy.  Some  of  these  structures 
are  shown  in  Fig.  2.  Many  of  these  also  have  a  high 
symmetry,  but  different  from  any  simple  lattice  sym¬ 
metry.  Generally,  up  to  the  size  of  12  atom  cluster,  all 
the  atoms  are  evenly  distributed  on  a  shell  having  roughly 
the  same  distance  from  the  center  of  the  cluster  (there  is 
no  atom  in  the  center  of  the  cluster).  The  13  atom  cluster 
is  an  icosahedron  with  an  atom  in  the  center,  but  again  in 
the  14  atom  cluster  all  the  atoms  are  on  one  shell,  the 
center  being  empty.  From  15  to  21  atom  clusters  there  is 
one  atom  in  the  center  and  the  rest  of  the  atoms  are 
evenly  distributed  on  the  surface,  nearly  the  same  dis¬ 
tance  from  the  center.  Larger  clusters  start  to  have  more 
atoms  also  inside  the  cluster,  but  the  surface  always 
seems  to  take  as  spherical  a  shape  as  possible. 

From  Figure  1  it  is  obvious  that  for  only  a  few  cluster 
sizes  does  the  lattice  structure  come  even  close  to  the 
binding  energy  of  the  fully  relaxed  optimal  configuration. 
The  large  variation  as  a  function  of  the  cluster  size  disap¬ 
pears  almost  totally  when  the  ions  are  allowed  to  relax 
from  the  lattice  sites  to  the  optimal  ground  state 
configuration.  The  remaining  variation  is  slightly  smaller 
than  the  variation  of  the  purely  electronic  part  of  the 
energy  in  the  jellium  model17  (coming  essentially  from 
the  first  three  terms  of  Eq.  (2)).  This  is  clearly  seen  in 
Fig.  3  where  the  second  derivative  of  the  energy  versus 
cluster  size  curve  is  shown  for  the  Madelung  energy  and 
for  the  energy  of  the  jellium  model.  This  result  has  the 
important  consequence  that  even  if  the  actual  interactions 
between  the  electrons  and  the  ions  are  taken  into  account, 
magic  numbers  may  still  be  determined  by  the  filling  of 
the  single  particle  energy  levels  in  a  spherical  potential. 
The  results  of  surface  energy  calculations  in  a  related 
model  show  that  the  structure  dependence  will  be  reduced 
if  the  electron  density  is  allowed  to  relax  from  the  jellium 
result14.  It  is  then  expected  that  the  variation  of  the 
structural  part  of  the  energy  is  overestimated  if  only  the 
Madelung  energy  is  considered.  The  noncrystalline  clus¬ 
ter  structures  seem  to  persist  up  to  cluster  sizes  of  several 
hundreds  of  atoms.  This  is  indicated  by  the  results  of 
Ihigues  et  alx  which  show  that  the  binding  energies  of  the 
crystalline  structures  remain  lower  than  those  of  the 
optimal  structures  even  at  cluster  sizes  of  more  than  90 


atoms.  Also  the  calculations  of  the  reconstruction  of  sim¬ 
ple  metal  surfaces18  indicate  that  the  registration  shift 
extends  several  layers  below  the  surface  and  originates 
from  the  same  electrostatic  forces  as  the  noncrystalline 
cluster  structures.  (The  cluster  size  has  to  be  larger  than 
about  500  atoms  before  more  than  half  of  the  atoms  are 
not  in  the  surface  layer).  For  very  large  clusters,  of 
course,  the  strictly  spherical  shape  will  be  modified  by 
faceting  which  will  eventually  make  the  crystalline  struc¬ 
ture  more  stable  than  the  (disordered)  spherical  one. 

In  conclusion,  we  have  shown  that  the  main  part  of 
the  structural  dependence  of  the  total  energy  of  small 
simple  metal  clusters  comes  from  the  classical  Madelung 
energy.  The  calculations  for  ground  state  structures  show 
that  this  Madelung  energy  varies  as  a  function  of  the 
cluster  size  less  than  the  electronic  contribution  which 
shows  the  shell  structure.  This  result  gives  further  under¬ 
standing  why  the  magic  numbers  can  be  understood  in 
terms  of  a  spherical  jellium  model,  and  predicts  that  clus¬ 
ter  structures  differ  from  crystalline  structures  up  to  clus¬ 
ter  sizes  of  hundreds  of  atoms. 

Acknowledgements  -  The  author  would  like  to  thank  N.W. 
Ashcroft,  J.W.  Wilkins  and  J.H.  Davies  for  useful  discus¬ 
sions.  This  work  was  supported  by  the  Office  of  Naval 
Research. 


REFERENCES 


1.  W.D.  Knight.  K.  Clemenger,  W.A.  de  Hcer,  W.A. 
Saunders.  MY.  Chow  and  ML.  Cohen.  Phvs.  Rev. 
Lett.  52,2141  (1984). 

2.  W.D.  Knight,  W.A.  de  Heer,  K.  Clemenger  and 
W.A.  Saunders,  Solid  State  Commun.  53,  445 
(1985). 

3.  M.Y.  Chow,  A.  Cleland  and  ML.  Cohen,  Solid  State 
Commun.  52,  645  (1984). 

4.  M.  Manninen,  R.M.  Nieminen  and  M.J.  Puska,  Phvs. 
Rev.  B  (in  press). 

5.  W.D.  Knight,  K.  Clemenger,  W.A.  de  Haar,  and 
W.A.  Saunders,  Phys.  Rev.  B  31,  2539  (1985). 

7.  J.L.  Martins,  J.  Buttet,  and  R.  Car,  Phys.  Rev.  B  31, 
1804  (1985). 

8.  M.P.  Iniguez.  J.A.  Alonso,  and  L.C.  Baibas,  Solid 
State  Commun.  57,  85  (1986). 

9.  J.L.  Martins,  R.  Car,  and  J.  Buttet,  Surf.  Sci.  106, 
265  (1981). 

10.  N.D.  Lang  and  W  Kohn,  Phys.  Rev.  B  1,  4555 
(1970). 

11.  I.V.  Abarenkov  and  V.  Heine,  Phil.  Mag.  12,  529 
(1965). 

12.  N.W.  Ashcroft,  Phys.  Lett.  23 A,  48  (1966). 

13.  The  electronic  structure  is  calculated  in  a  spherically 
averaged  potential  which  depends  on  the  sites  of  the 
pseudopotentials.  Due  to  the  averaging  the  structure 
dependence  is  diminished. 

14.  R.  Monnier  and  J.P.  Perdew,  Phys.  Rev.  B  17,  2595 
(1978). 

15.  R.N.  Barnett,  U.  Landman,  and  C.L.  Cleveland, 
Phys.  Rev.  Lett.  51,  1359  (1983). 

16.  M.  Manninen  and  R.M.  Nieminen,  J.  Phys.  F  8, 
2243  (1978). 

17.  B.K.  Rao,  P.  Jena,  and  M.  Manninen,  Phys.  Rev.  B 
32,477  (1985). 


7 


FIGURE  CAPTIONS 


Fig.  1.  Calculated  Madelung  energies  as  a  function  of  the 
cluster  size.  The  upper  curve  shows  the  results  for  the 
fully  relaxed  structures  and  the  lower  curve  is  the  result 
for  the  most  stable  bcc,  fee  or  hep  structures.  The 
Madelung  energy  for  the  bulk  bcc  and  fee  lattices  is  also 
shown  in  the  figure.  ESi  is  given  in  atomic  units  for  a 
density,  n0,  corresponding  to  r5=l. 


Fig.  2.  Examples  of  cluster  structures  determined  by 
minimizing  the  Madelung  energy. 


Fig.  3.  The  variation  of  the  Madelung  energy  (in  atomic 
units)  as  a  function  of  the  cluster  size  (a)  as  compared  to 
the  variation  of  the  total  energy  of  the  jellium  model17 
(b).  The  results  are  for  lithium  (r5=3.25). 
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Table  1.  Optimal  structures  of  spherical  clusters  assum¬ 
ing  them  to  have  bcc,  fee,  or  hep  structures.  The  center 
of  the  cluster  has  been  chosen  to  be  lattice  site  (lat),  mid¬ 
point  of  two  atoms  (mid),  tetrahedral  interstitial  site  (ter), 
octahedral  interstitial  site  ( oct ),  or  the  center  of  an  equi¬ 
lateral  triangle  (tri). 


Size 

Lattice(center) 

1 

2 

bcc(mid) 

3 

fcc(tri) 

4 

fcc(tet) 

5 

fcc(oct) 

6 

fcc(oct) 

7 

fcc(oct) 

8 

fcc(oct) 

9 

fcc(oc0 

10 

hcp(mid) 

11 

bcc(lat) 

12 

fcc(lat) 

13 

fcc(lat) 

14 

bcc(lat) 

15 

bcc(lat) 

16 

bcc(lat) 

17 

hcp(lat) 

IS 

bcc(lat) 

19 

fcc(lat) 

20 

hcp(tet) 

Size 

Lattice(center) 

21 

hcp(tet) 

22 

hcp(tet) 

23 

hep^mid) 

24 

hcp(tet) 

25 

hcp(tri) 

26 

hcp(tri) 

27 

hcp(tri) 

2S 

hcp(tet) 

29 

hcp(tet) 

30 

hcp(tet) 

31 

bcc(tet) 

32 

hcp(tet) 

33 

hcp(oct) 

34 

bcc(mid) 

35 

hcp(oct) 

36 

fcc(oct) 

37 

fcc(oct) 

3S 

t'cc(oct) 

39 

fcc(oct) 

40 

fcei,oct) 

